Coupling/decoupling between translational and rotational dynamics 
in a supercooled molecular liquid 
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We use molecular dynamics computer simulations to investigate the coupling/ decoupling between 
translational and rotational dynamics in a glass-forming liquid of dumbbells. This is done via 
a careful analysis of the a-relaxation time Tq* of the incoherent center-of-mass density correlator 
at the structure factor peak, the a-relaxation time T2 of the reorientational correlator, and the 
translational (D t ) and rotational (D r ) diffusion constants. We find that the coupling between the 
relaxation times t5 and T2 increases with decreasing temperature T, whereas the coupling decreases 
between the diffusivities Dt and D r . In addition, the T-dependence of Dt decouples from that of 
1/t"2, which is consistent with previous experiments and has been interpreted as a signature of the 
"translation-rotation decoupling." We trace back these apparently contradicting observations to the 
dynamical heterogeneities in the system. We show that the decreasing coupling in the diffusivities 
Dt and D r is only apparent due to the inadequacy of the concept of the rotational diffusion constant 
for describing the reorientational dynamics in the supercooled state. We also argue that the coupling 
between Tq* and T2 and the decoupling between Dt and I/V2, both of which strengthen upon cooling, 
can be consistently understood in terms of the growing dynamic length scale. 

PACS numbers: 64.70.pm, 61.20.Ja, 61.20.Lc, 61.25.Em 



One of the puzzling features of glass- forming systems is 
that upon cooling the translational dynamics appears to 
decouple from the rotational dynamics [J 0, S H, @ • At 
high temperatures, the translational and rotational dif- 
fusion coefficients, D t and D r , respectively, are inversely 
proportional to rj/T, where 77 is the viscosity and T the 
temperature, in agreement with the Stokes-Einstein (SE) 
and the Stokes-Einstein-Debye (SED) relation. Below 
approximately 1.2 T g , with T g the glass transition tem- 
perature, the SE relation is replaced by a fractional rela- 
tion D t oc (ry/T) - * with £ < 1, implying a more enhanced 
translational diffusion than predicted by the viscosity, 
whereas D r remains proportional to (rj/T)^ 1 down to 
T g [ll, 0). On the other hand, recent computer simu- 
lations give evidence for a stronger correlation between 
translational and rotational mobilities at lower T 0, 0] , 
which is apparently in contradiction with the mentioned 
decoupling between Dt and D r . In addition, an enhance- 
ment of rotational diffusion relative to translation upon 
cooling is reported in these simulations, which is a trend 
opposite to that observed in experiments [H, 0] . In this 
Letter, we report computer simulation results for gain- 
ing more insight and unified understanding of these cou- 
pling/decoupling phenomena in the supercooled state. 

The system we consider is the binary mixture of rigid, 
symmetric dumbbell molecules, denoted as AA and BB 
dumbbells, studied in Ref. Q. Each molecule consists of 
two identical fused Lennard- Jones (L J) particles of type 
A or B having the same mass m, and their bond lengths 
are denoted by Iaa and Ibb- The interaction between 
two molecules is given by the sum of the LJ interactions 
V a fj{r) between the four constituent sites, with the LJ 



parameters e a p and o nii for a, (3 € {A,B} taken from 
Ref. [§], which are slightly modified so that V a p(r) and 
V^Jr) are zero at the cutoff r cut = 2.5a a p (see Ref. [1] 
for details). Bond lengths are specified by a parameter 
£ = Iaa/oaa = Ibb/obb, and a sufficiently large value 
C = 0.8 is chosen so that anomalous reorientational dy- 
namics caused by the so-called type-A transition is ab- 
sent |8j, [lO(. The number of AA and BB dumbbells is 
N AA = 800 and N B b = 200. In the following all quanti- 
ties are expressed in reduced units with the unit of length 
ctaa, the unit of energy 6aa (setting ks = 1), and the 
unit of time {ma\ A / tAA) 1 ^ 2 ■ Standard molecular d 
namics simulations have been performed as in Ref. 
with the cubic box of length L = 10.564 for 2.0 < T < 10. 
The longest runs were 2 x 10 9 time steps, and we per- 
formed 16 independent runs to improve the statistics. 
Such long simulations were necessary to reach below the 
critical temperature T c of the mode-coupling theory [llj , 
which we estimate as T c w 2.10 based on an analysis 
similar to the one done in Ref. Q. In experiments T c 
is found to be « 1.2 T g [l2|, i.e., close to the tempera- 
ture at which the aforementioned decoupling sets in [l[. 
We also introduce the onset temperature, found to be 
Tonsot ~ 4.0, below which correlators exhibit the two-step 
relaxation, a characteristic feature of the glassy dynamics 
in which molecules are temporarily caged by their neigh- 
bors. Hereafter, all quantities used in our discussion refer 
to those for A A dumbbells, and the subscript A A shall 
be dropped for notational simplicity. 

In the present study we will characterize the trans- 
lational dynamics by the incoherent center-of-mass den- 



sity correlator F°(t) = (l/N)J2j< 
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FIG. 1: Correlation coefficient between rHj and T2j as a 
function of 1/T. The Inset is a scatter plot of T2j versus t£*j 
for T = 2.0 on double logarithmic scales. 



the rotational dynamics by Ci(t) = (1/N) Y^, j(PeWj{t) 
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FIG. 2: (Color online) Log-log plot of the translational dif- 
fusion coefficient Dt (circles, right scale), the rotational dif- 
fusion coefficient D r (diamonds, right scale, shifted so that 
Dt and D r agree at high T), and the inverse of the rotational 
relaxation time 1/t2 (filled squares, left scale) versus the a- 



% (°)]) ■ Here ?f (*) and % (<) denote the center-of-mass re i axat ion time r 9 c . . The dotted straight lines refer to (re- 



position and the orientation of the jth molecule at time t, 
respectively, and Pi is the Legendre polynomial of order 
I. The a-relaxation times and ti shall be defined via 
if Of) = 0.1 and C^n) = 0.1. 

Previous simulation studies have used various ways 

to classify particle mobility @, S El Q (H. in this 

Letter, translational and rotational mobilities of individ- 
ual molecule shall be classified in terms of the first pas- 
sage times, Tqj and Tg t j, at which individual- molecule 

quantities e % ^ r i ^~ r i (averaged over q having the 
same modulus q) and Pi [ej (t) ■ ej (0)] become zero for the 
first time. We confirmed that oc (1/N) ^ • j and 
ti oc (1/N) J2j T e,j hold, and hence, r^j and tij can also 
be regarded as individual molecule's relaxation times. In 
the following, we will mainly refer to the structure fac- 
tor peak position q* = 8.0 and t = 2 because of their 
experimental significance [j], @, [l6[ . 

The Inset of Fig. Q] shows a scatter plot of raj versus 
r^ ■ for T = 2.0. We recognize that there is a strong 
correlation between these two quantities, in particular 
for the molecules that move/rotate quickly. From such 
a scatter plot we can calculate the correlation coefficient 
and its T-dependence is shown in the main panel of Fig. [1] 
We see that the correlation between translation and ro- 
tation increases quickly with decreasing T, showing the 
strong coupling between these two types of motions at 
low temperatures, in agreement with the findings in re- 
cent simulation studies (y, l7( ■ 

We also mention that the molecules that have a high 
translational mobility are clustered in space, in agree- 
ment with the findings for other glass formers 0, H llo| , 
and the same is true for those rotating quickly. This 
result is in accord with the correlation shown in Fig. [TJ 

We next consider the translational (D t ) and rotational 



with the exponents x cited in the figure. The arrows indicate 
the locations of T onsct and T c . Inset: D t (circles) versus r^, 
from the main panel is compared with (I/t^L j) (solid line) and 
1/Tq for q = 1.2 (dashed line), which are vertically shifted so 
that they agree with D t at high T. 



(D r ) diffusion coefficients. D t is determined via D t — 
limt^oo Ar^ (t) /6t of the center-of-mass mean-squared 
displacement Ar 2 c (t) = (1/N) XV([ff(t) - r?(0)] 2 ). D r 
is calculated using the Einstein formulation [17] . i.e., via 
D r = lim^oo A(j> 2 (t) / At of the mean-squared angular 
displacement 
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Here A<^ )Q (i) — L dt' cjj tCt (t') with u>j tCt (t) denoting the 
a-component of the angular velocity of the jth molecule. 
Note that for linear molecules, the summation over a is 
for the two directions X and Y perpendicular to the Z 
axis chosen along £?j(0) [Hj]. 

Figure O shows the log- log plot of D t (circles) and D r 
(diamonds) versus r^ . Thus, this graph mimics Fig. 5 of 
Ref. [l[ where D t and D r for orthoterphenyl are traced 
as a function of rj/T. We have used substitute 
for rj/T since accurate calculation of r\ is computation- 
ally challenging. However, this approximation should be 
quite accurate since the a-relaxation time at the struc- 
ture factor peak position is known to track rj/T (l6| . 

From Fig. [2] we see that D t oc (r^) _1 holds at high 
T, whereas it is replaced by a fractional relation D t oc 
( T C)-o.83 for t < T c . This anomalous behavior (the 
breakdown of the SE relation) and its occurrence near 
T c w 1.2 T g are in agreement with experiments [lj]. For 
D r we recognize that D r oc (r^) _1 , and hence D r oc D t , 
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FIG. 3: (Color online) (a) The reorientational correlator 
Ci(t) (solid lines) and the one calculated within the Gaussian 
approximation C§{t) (dashed lines) versus log 10 t for T = 5.0, 
3.0, 2.39, and 2.18 (from left to right), (b) The average orien- 
tation 9(t) (left scale) and the logarithm of the mean-squared 
angular displacement log 10 A(f> 2 (t) (right scale) for T = 2.18. 

hold at high T, but that this quantity strongly decouples 
from D t for T < T onsct . That the T-dependence of D r 
is weaker than the one of D t is in contradiction to the 
trend observed in experiments [lj], but is in agreement 
with other recent simulation studies @, 0] • 

Let us notice that the experiment in Ref. [l[ does not 
directly probe D r , but determines D r from the measure- 
ment of T2 via the relation D r cx 1/t 2 predicted by the 
Debye model. Thus, what is reported as D r in Ref. [1] 
corresponds in fact to l/r 2 . In Fig. [2] we have also in- 
cluded our simulation results for 1/t2 and we see that 
T2 correlates well with , in accord with our previous 
discussion for the relaxation times, although a slight de- 
viation from the strict proportionality is discernible for 
T < T c . Thus, by identifying rfji with r//T, our simula- 
tion result for 1 /n is in accord with the experiment [H in 
which one finds 1/t2 oc (rj/T)^ 1 down to T g . Therefore, 
the decoupling of D r from I/V2 shown in Fig. [2] only indi- 
cates that the Debye model breaks down for T < T onset , 
and hence, I/V2 should not be considered as being pro- 
portional to D r in the supercooled state. 

To understand why D t and D r behave so differently, 
we analyze the breakdown of the Debye model in more 
detail. Within the Gaussian approximation one has 
Cf(t) = expH(£+ l)A(f> 2 {t)/4} [18]. The long-time be- 
havior is given by Cf (t) = exp[— £(£ + l)D r t], which 
agrees with the prediction from the Debye model. From 
Fig. [5] one therefore expects that the Gaussian approx- 
imation breaks down for T < T onsct « 4.0, and this is 
indeed the case as shown in Fig. [3K 

Figure [3^ also implies that for T < T onsct there ex- 



ists a time window during which A(j> 2 (t) increases sig- 
nificantly although the orientation of a molecule hardly 
changes. The existence of such a window is demon- 
strated in Fig. [3Jd for a representative temperature, from 
which one recognizes that A(f> 2 (t) starts to increase no- 
ticeably at log 10 1 k, 1 whereas the average orientation 
9{t) = {l/N) J2j ( cos _1 [<5j (t) ■ ej(0))) remains nearly con- 
stant up to log 10 1 « 2. 

The appearance of the plateau in C2(t) and 9(t) for 
T < T onset reflects the librational motion of a molecule 
trapped inside the cage, with the orientation ej(t) of 
the molecule remaining close to ej(0). (For the transla- 
tions! degrees of freedom this corresponds to rattling of 
the center-of-mass position.) If the average orientation 
during the librational motion coincides with e^fO), this 
gives rise to a diffusive movement of A.(f>j t z(t) [17] . which 
however does not contribute to A(f> 2 (t), see Eq. JT]). If, 
on the other hand, the average orientation is tilted from 
e*j(0), one can show that A.(pj t x(t) and A0j.y(£) also ex- 
hibit a diffusive movement. Thus, a diffusive increase 
in A(j) 2 (t) occurs even if the molecule has not yet reori- 
ented and is still performing the librational motion inside 
the cage, and this effect explains the difference between 
log 10 A(j) 2 (t) and 0(t) at intermediate times as seen in 
Fig. [3Jd. D r from the Einstein formulation is thus more 
affected by the librational motion than the real "reorien- 
tational" dynamics, leading to a spurious "decoupling" 
between D t and D r . Since there is no other proper way 
to define D r , we conclude that the concept of the rota- 
tional diffusion constant is inadequate for describing the 
reorientational dynamics in the supercooled state. 

Now let us turn our attention to the decoupling be- 
tween 1/Tq, and D t , and between 1/t 2 and D t , which oc- 
curs near and below T c (see Fig. [5]). The latter is known 
as the "translation-rotation decoupling." It was conjec- 
tured that such a decoupling arises because 1/t 2 and D t 
reflect different moments of the distribution of relaxation 
times, i.e., l/r 2 oc l/(r 2j ) while D t cx (1/t^-) >]. The 
essential ingredient of this conjecture is the assumption 
that local translational mobility is proportional to local 
rotational mobility 4], which is justified from our simula- 
tion (see Fig.Q]). However, while T2 cx (t 2j ) has been con- 
firmed from our simulation, the T-dependence of (l /r^ ■) 
is not in accord with that of D t as demonstrated in the 
Inset of Fig. [2l On the other hand, as also shown there, 
the inverse of the a-relaxation time for q = 1.2 is 

found to track D t in the simulated T range. 

That D t oc (r g c =1 while D t oc (r^ )-° 83 for T < T c 
can be understood in terms of t he g rowing onset length 
scale ^onsct of Fickian diffusion [ijj [2(|. Here, /onset is 
defined so that the relaxation time for 27r/g < l onsc t 
decouples from the diffusivity D t , and can be estimated 
from the ratio B% of the product q 2 D t r^ to the one at 
some reference temperature T rc f with a procedure de- 
tailed in Ref. [2C|. (In the present work we have cho- 
sen T rc f = 5.0 and a criterion R™ = 1.2 for determin- 
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FIG. 4: (Color online) Correlation coefficients between t£* j 
and Tg^'s (solid lines) for q = 4.2, 2.4, and 1.2 versus 1/T. 
The dashed line denotes the correlation coefficient between 
T2j and riLi.2 j. The Inset exhibits the length scale ionsot 
introduced in the text as a function of 1/T. 



ing Zonsct-) The resulting Z onS ct as a function of 1/T is 
presented in the Inset of Fig. 01 from which one finds 
Zonset « 3 for 1/T « 1/T C = 0.48. Thus, D t cx (r^^)- 1 
for T < T c since the associated length scale 27r/1.2 ps 
5.2 exceeds Z nset, whereas decouples from D t since 
2ir/q* ~ 0.8 is smaller than i nset- 

To what extent the dynamics occurring on different 
length scales are correlated is examined in the main panel 
of Fig. [H where correlation coefficients between rS • and 
Tqj's for q — 4.2, 2.4, and 1.2 are plotted with solid 
lines. The mentioned decoupling between and t^ =1 2 
is reflected in the small value (< 0.1) of the correspond- 
ing correlation coefficient. On the other hand, one in- 
fers from Fig. |4] that a significant correlation, character- 
ized, e.g., by the correlation coefficient exceeding 0.2, de- 
velops between rSij and t^L 4 2 j for 1/T > 0.4, where 
^onsct becomes larger than the associated length scale 
27r/4.2 ps 1.5. Similarly, one finds that the correlation 
coefficient between t^, • and r£L 2 4 j exceeds 0.2 only for 
1/T > 0.5, where Z onsot ' is larger than 2tt/2.4 ps 2.6. Thus, 
though originally introduced as the onset length scale of 
Fickian diffusion, Z onso t also provides a coherent length 
in the sense that the dynamics occurring on length scales 
smaller than Z onsc t develop significant correlations. 

Let us estimate the length scale associated with t 2 uti- 
lizing such a "coherent length" l onse t. One finds from 
Fig. Q] that the correlation coefficient between r£i a and 
T2.j exceeds 0.2 for 1/T > 0.4. According to the men- 
tioned interpretation of Zonsct, the dynamics responsible 
for Tq, and t 2 become correlated if both of their charac- 
teristic length scales become smaller than ionsct- Thus, 
the length scale associated with t 2 can be estimated as 
ps 1.5 which is determined by Z nsct at 1/T pa 0.4 (see 
Fig. [4]) . This explains why the correlation coefficient be- 
tween T2,j and T^L 12 ,j i s small 0.1) as shown with 



the dashed line in Fig. [4] i.e., why 1/t 2 decouples from 

Thus, the translation-rotation decoupling between D t 
and T2 should more properly be understood as the de- 
coupling of the dynamics occurring on different length 
scales, which arises due to the growing dynamic length 
scale. This notion also explains the decoupling between 
t^, and D t (the breakdown of the SE relation), as well 
as the translation-rotation "coupling" between and 
t 2 occurring on comparable length scales. 
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